AFRL-IF-RS-TR-2006-295 
Final  Technical  Report 
September  2006 


VERY  LARGE-SCALE  MULTIUSER  DETECTION 
(VLSMUD) 


Princeton  University 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


STINFO  FINAL  REPORT 


AIR  FORCE  RESEARCH  LABORATORY 
INFORMATION  DIRECTORATE 
ROME  RESEARCH  SITE 
ROME,  NEW  YORK 


NOTICE  AND  SIGNATURE  PAGE 


Using  Government  drawings,  speeifieations,  or  other  data  ineluded  in  this  doeument  for 
any  purpose  other  than  Government  proeurement  does  not  in  any  way  obligate  the  U.S. 
Government.  The  faet  that  the  Government  formulated  or  supplied  the  drawings, 
speeifieations,  or  other  data  does  not  lieense  the  holder  or  any  other  person  or  eorporation; 
or  eonvey  any  rights  or  permission  to  manufaeture,  use,  or  sell  any  patented  invention  that 
may  relate  to  them. 

This  report  was  cleared  for  public  release  by  the  Air  Force  Research  Laboratory  Rome 
Research  Site  Public  Affairs  Office  and  is  available  to  the  general  public,  including  foreign 
nationals.  Copies  may  be  obtained  from  the  Defense  Technical  Information  Center  (DTIC) 
(http://www.dtic.mil). 


AFRL-IF-RS-TR-2006-295  HAS  BEEN  REVIEWED  AND  IS  APPROVED  EOR 
PUBEICATION  IN  ACCORDANCE  WITH  ASSIGNED  DISTRIBUTION 
STATEMENT. 


EOR  THE  DIRECTOR: 


/s/ 


/s/ 


DAVID  HENCH 
Work  Unit  Manager 


WARREN  H.  DEBANY,  Jr. 

Technical  Advisor,  Information  Grid  Division 
Information  Directorate 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange,  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 

OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 

of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Service,  Directorate  for  Information  Operations  and  Reports, 

1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget, 

Papen/vork  Reduction  Project  (0704-0188)  Washington,  DC  20503. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

SEP  2006 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  TYPE 


3.  DATES  COVERED  (From  -  To) 

Apr  05  -  Aug  06 


5a.  CONTRACT  NUMBER 


VERY  LARGE  SCALE  MULTI-USER  DETECTION  (VLMUD) 


5b.  GRANT  NUMBER 

FA8750-05-2-0I92 


5c.  PROGRAM  ELEMENT  NUMBER 

35885G 


5d.  PROJECT  NUMBER 


VLSM 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Princeton  University 

Office  of  Research  and  Project  Administration 
New  South  Building,  Princeton  NJ  08544-0036 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


AFRL/IFGC 
525  Brooks  Rd 
Rome  NY  13441-4505 


12.  DISTRIBUTION  AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED.  PA#  06-< 


11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

AFRL-IF-RS-TR-2006-295 


14.  ABSTRACT 


This  document  is  the  final  report  for  award  FA8750-05-2-0192.  It  summarizes  findings  in  our  study  of  Very  Large  Scale  Multi-user 
Detection.  This  work  has  involved  two  distinct  efforts:  the  first  addressed  the  analysis  of  iterative  channel  estimation  and  multi-user 
detection  in  multi-path  DS-CDMA  channels,  and  the  second  addressed  the  issue  of  energy  efficiency  in  multi-hop  CDMA  networks. 


15.  SUBJECT  TERMS 

Iterative  channel  estimation,  multi-path  DS-CDMA  channels,  multi-hop  CDMA  networks,  multi-user  detection 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 

OF  PAGES 

ISa.  NAME  OF  RESPONSIBLE  PERSON 

David  Flench 

a.  REPORT 

U 

b.  ABSTRACT 

U 

c.  THIS  PAGE 

U 

UL 

45 

ISb.  TELEPHONE  NUMBER  [Include  area  code) 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


TABLE  OF  CONTENTS 


1 .  Summary . 1 

2.  Introduction . 1 

3.  Methods,  assumptions  and  procedures . 1 

4.  Results  and  discussion . 2 

5.  Conclusions . 2 


Appendix  A:  Performance  Analysis  of  Iterative  Channel  Estimation  and  Multiuser  Detection 

in  Multipath  DS-CDMA  Channels . 3 

Appendix  B;  Energy  Efficiency  in  Multi-hop  CDMA  Networks:  A  Game  Theoretic  Analysis* 
. 36 


i 


1.  SUMMARY 

Our  work  in  this  area  has  involved  two  distinct  efforts:  the  first  addresses  the  analysis  of  iterative 
channel  estimation  and  multiuser  detection  in  multipath  DS-CDMA  channels  [1],  and  the  second 
addresses  the  issue  of  energy  efficiency  in  multi-hop  CDMA  networks  [2], 

2.  INTRODUCTION 

Multiuser  detection  refers  to  the  data  demodulation  in  multiple-access  communication  networks. 
It  is  a  key  technology  in  the  development  of  high-capacity  wireless  communications  and  sensor 
networks.  This  study  was  concerned  with  multiuser  detection  in  networks  with  very  large 
numbers  of  nodes,  i.e.,  very  large  scale  multiuser  detection.  Our  work  has  focused  on  two  aspects 
of  this  problem.  One  of  these  is  in  establishing  effective  methods  of  analysis  for  multiuser 
receivers  that  iterate  between  multiuser  detection  (and  channel  decoding)  and  channel  estimation 
[1].  Channel  estimation  is  necessary  in  wireless  receivers  because  of  the  dynamic  fading  that 
characterizes  wireless  communication  channels.  Although  such  iterative  receivers  have  been 
proposed  before,  their  analysis  (and  thus  their  effectiveness)  has  been  an  open  question.  The 
second  aspect  of  this  problem  is  that  of  the  effects  of  multiuser  detection  on  the  energy  efficiency 
of  multi-hop  networks  [2],  which  are  of  considerable  interest  in  a  variety  of  applications  of 
interest  to  the  Air  Force.  Energy  efficiency  is  often  a  primary  consideration  in  design  and 
deployment,  and  multiuser  detection  can  have  considerable  effects  on  this  efficiency.  The 
purpose  of  this  work  is  to  develop  methods  that  can  be  used  to  characterize  these  effects  in 
support  of  design  criteria  for  the  such  networks. 

3.  METHODS,  ASSUMPTIONS  AND  PROCEDURES 

This  work  is  described  in  detail  in  the  two  papers  [1,2],  copies  of  which  are  attached  as 
Appendices.  Generally  speaking  the  primary  assumption  is  that  of  a  direct-sequence  code¬ 
division  multiple-access  (DS-CDMA)  network  with  a  very  large  number  of  nodes.  The  principal 
analytical  method  used  to  examine  the  very  large-scale  case,  is  the  large-system  limit  of 
DSCDMA  systems,  in  which  the  spreading  gain  and  the  number  of  nodes  both  increase  without 
bound  while  their  ratio  (the  “system  load”)  remains  constant.  This  analytical  device  allows  for 
closed-form  analysis  in  both  of  the  problems  addressed.  To  analyze  the  iterative  channel 
estimation  and  multiuser  detection  receiver  in  [1]  the  main  additional  techniques  are  those  used 
to  study  the  convergence  of  iterative  algorithm.  To  study  the  energy  efficiency  of  multi-hop 
networks  in  [2]  a  game-theoretic  approach  is  used,  in  which  the  terminals  in  the  network  work 
are  view  as  “economic”  agents  competing  for  radio  resources  to  transmit  their  messages  as 
efficiently  as  possible  (i.e.,,  with  maximal  bits  transmitted  successfully  per  joule  of  battery 
energy). 
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4,  RESULTS  AND  DISCUSSION 

The  results  are  deseribed  in  detail  in  the  aeeompanying  publieations  [1,2].  In  the  problem  of 
iterative  ehannel  estimation  and  multiuser  deteetion  of  [1],  the  results  are  in  the  form  of 
eharaeterization  of  the  eonditions  (in  terms  of  system  load  and  SNR)  under  whieh  the  iterative 
methods  eonverge,  and  the  resulting  performanee.  In  the  problem  of  energy  efficiency  in  multi¬ 
hop  networks  of  [2],  the  results  are  in  the  form  of  the  energy  efficiency  per  node  (bits-per-joule) 
as  a  function  of  the  type  of  multiuser  detector  used  by  the  terminals. 

5.  CONCLUSIONS 

For  detailed  conclusions,  please  see  the  accompanying  publications,  [1,2],  in  the  appendix. 
Briefly,  it  is  seen  in  [1]  that  iteration  can  substantially  improve  upon  the  traditional 
(noniterative)  methods  of  multiuser  receiver  implementation.  In  [2],  it  is  seen  that  the  use  of 
multiuser  detection  can  be  a  major  determinant  in  the  energy  efficiency  of  multi-hop  multiuser 
systems. 

APPENDIXES 

[1]  H.  Li,  S.  Betz  and  H.  V.  Poor,  “Performance  Analysis  of  Iterative  Channel  Estimation  and  Multiuser 
Detection  in  Multipath  DS-CDMA  Channels,”  IEEE  Transactions  on  Signal  Processing,  accepted  for 
publication.  [Originally  submitted  December  14,  2005;  revision  submitted  March  13,  2006;  accepted 
August  15,  2006.] 

[2]  S.  Betz  and  H.  V.  Poor,  “Energy  Efficiency  in  Multi-hop  CDMA  Networks:  A  Game  Theoretic 
Analysis,”  (invited  paper).  Proceedings  of  the  Workshop  on  Multi-Layer  Modelling  and  Design  of  Multi- 
Hop  Wireless  Networks  (MEMO ’06),  Minneapolis,  MN,  July  12  -  15,  2006. 
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Appendix  A 

Performance  Analysis  of  Iterative  Channel 
Estimation  and  Multiuser  Detection  in  Multipath 

DS-CDMA  Channels 

Husheng  Li,  Sharon  M.  Betz  and  H.  Vincent  Poor 


Abstract 

This  paper  examines  the  performance  of  decision  feedback  based  iterative  channel  estimation  and  multiuser  de¬ 
tection  in  channel  coded  aperiodic  DS-CDMA  systems  operating  over  multipath  fading  channels.  First,  explicit  ex¬ 
pressions  describing  the  performance  of  channel  estimation  and  parallel  interference  cancellation  based  multiuser 
detection  are  developed.  These  results  are  then  combined  to  characterize  the  evolution  of  the  performance  of  a  system 
that  iterates  among  channel  estimation,  multiuser  detection  and  channel  decoding.  Sufficient  conditions  for  conver¬ 
gence  of  this  system  to  a  unique  fixed  point  are  developed. 


I.  Introduction 

Direct  sequence  code  division  multiple-access  (DS-CDMA)  has  been  selected  as  the  fundamental  sig¬ 
naling  technique  for  third  generation  (3G)  wireless  communication  systems,  due  to  its  advantages  of  soft 
user  capacity  limit  and  inherent  frequency  diversity.  However,  it  suffers  from  multiple-access  interference 
(MAI)  caused  by  the  non-orthogonality  of  spreading  codes,  particularly  for  heavily  loaded  systems.  There¬ 
fore,  techniques  for  mitigating  the  MAI,  namely  multiuser  detection,  have  been  the  subject  of  an  intensive 
research  effort  over  the  past  two  decades.  It  is  well  known  that  multiuser  detection  can  substantially  sup¬ 
press  MAI,  thus  improving  system  performance.  Maximum  likelihood  (ML)  multiuser  detection  [28]  was 
proposed  in  the  early  1980s,  and  achieves  the  optimal  performance  at  the  cost  of  prohibitive  computational 

Husheng  Li  is  with  Qualcomm  Inc.,  San  Diego,  CA,  92121,  USA  (email:  hushengl@qualcomm.com).  Sharon  M.  Betz  is  with  Department  of 
Electrical  Engineering,  Princeton  University,  Princeton,  NJ  08544,  USA  (email:  sbetz@princeton.edu).  H.  Vincent  Poor  is  with  Department  of 
Electrical  Engineering,  Princeton  University,  Princeton,  NJ  08544,  USA  (email:  poor@princeton.edu).  This  research  was  supported  by  the  Air 
Force  Research  Laboratory  under  Cooperative  Agreement  No.  FA8750-05-2-01  92. 
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cost  when  the  number  of  users  is  large.  For  practical  implementation,  suboptimal  algorithms,  such  as  the 
linear  minimum  mean  square  error  (LMMSE)  detector  [21]  or  decorrelator  [29],  allow  a  tradeoff  between 
complexity  and  performance.  It  should  be  noted  that  the  teehnique  of  multiuser  detection  is  being  applied  in 
existing  CDMA  systems,  such  as  EV-DO  Revision  A  systems  [12]. 

In  recent  years,  the  turbo  principle,  namely  the  iterative  exehange  of  soft  information  among  different 
bloeks  in  a  eommunication  system  to  improve  the  system  performanee,  has  been  applied  to  combine  mul¬ 
tiuser  detection  with  channel  deeoding  [1] [22] [24] [26] [27] [31].  In  sueh  turbo  multiuser  deteetors,  the  out¬ 
puts  of  channel  decoders  are  fed  back  to  the  multiuser  detector,  thus  enhancing  the  performance  iteratively. 
Turbo  multiuser  deteetion  based  on  the  maximum  a  posteriori  probability  (MAP)  detection  and  decoding 
criterion  has  been  proposed  in  [30]  [31]  together  with  a  lower  eomplexity  teehnique  based  on  interferenee 
cancellation  and  LMMSE  filtering.  Eurther  simplifieation  is  obtained  by  applying  parallel  interferenee  can¬ 
cellation  (PIC)  [1]  for  multiuser  deteetion,  where  the  decisions  of  the  deeoders  are  directly  subtracted  from 
the  original  signal  to  eaneel  the  MAI. 

Practical  wireless  eommunieation  systems  usually  experience  fading  channels,  whose  state  information 
is  unknown  to  the  reeeiver.  Thus  practieal  systems  need  to  consider  detection  and  decoding  with  uncertain 
channel  state  information.  In  the  context  of  short  code  CDMA  systems,  blind  multiuser  deteetion  can  be 
accomplished  without  explicit  channel  estimation  by  using  subspace  and  other  teehniques  [32].  An  alter¬ 
native  reeeiver  strueture  adopts  an  explieit  channel  estimation  bloek  and  carries  out  the  decoding  with  the 
corresponding  ehannel  estimate.  In  systems  without  decision  feedback,  the  channel  estimation  block  is  cas¬ 
caded  with  the  decoder  and  operates  as  a  front  end  for  the  subsequent  blocks.  With  such  a  receiver  strueture, 
the  ehannel  estimates  ean  be  obtained  with  training  symbols  [6]  or  with  blind  estimation  algorithms  [33]. 
Explicit  expressions  for  the  performanee  of  such  channel  estimation  schemes  are  given  in  [17]  and  the  eor- 
responding  impact  on  multiuser  detection  is  discussed  in  the  large  system  limit  in  [9]  and  [18].  In  systems 
with  decision  feedback,  the  deeisions  of  the  decoder  are  fed  baek  to  the  ehannel  estimator  to  enhance  its 
performanee.  In  such  systems,  the  ehannel  estimator  and  the  decoder  ean  operate  either  simultaneously  [25] 
or  sueeessively  [7]  [13]  [23].  An  example  of  the  former  strategy  applied  to  ML  sequenee  deteetion  in  uncer¬ 
tain  environments  is  proposed  in  [25];  called  per-survivor  processing,  tentative  deeisions  are  immediately 
fed  baek  to  the  ehannel  estimation  algorithm  and  the  corresponding  estimates  are  used  for  the  detection 
of  future  symbols.  In  the  latter  strategy,  the  decisions  are  fed  baek  only  when  the  entire  eurrent  decoding 
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procedure  is  finished.  For  example,  in  [13],  an  expectation  maximization  (EM)  channel  estimation  algo¬ 
rithm,  combined  with  successive  interference  cancellation,  is  proposed.  Joint  channel  estimation  and  data 
detection  algorithms  for  uncoded  single- antenna  and  multiple-antenna  systems  are  discussed  in  [8]  and  [7], 
respectively.  In  channel  coded  systems,  iteration  can  achieve  better  performance  when  the  turbo  principle  is 
applied,  due  to  the  redundancy  introduced  by  the  code  structure.  In  [23],  an  iterative  algorithm  is  proposed 
and  analyzed  for  channel  estimation  and  decoding  of  low-density  parity-check  (LDPC)  coded  quadrature 
amplitude  modulation  (QAM)  systems. 

In  this  paper,  we  consider  channel-coded  CDMA  systems  operating  over  multipath  fading  channels  whose 
channel  state  information  is  unknown  to  the  receiver.  To  demodulate  and  deeode  such  systems,  we  apply 
the  turbo  principle  to  both  channel  estimation  and  multiuser  detection.  As  shown  in  Figure  1 ,  we  consider 
a  receiver  that  feeds  back  decisions  from  channel  decoders  to  both  an  ML  channel  estimator  and  a  PIC 
multiuser  detector.  The  iteration  is  initialized  with  training  symbol  based  channel  estimation  and  a  non¬ 
iterative  multiuser  detection.  The  reeeiver  structure  is  similar  to  those  proposed  in  [2]  [15]  [20].  However, 
this  paper  is  focused  mainly  on  the  performance  analysis  of  such  structures  using  semi-analytic  methods.  We 
analyze  the  contributions  to  the  variance  of  the  channel  estimation  error  due  to  noise  and  deeision  feedback 
error,  and  the  variance  of  the  residual  MAI  after  PIC.  We  then  use  this  analysis  to  describe  the  deeoding 
process  as  an  iterative  mapping.  We  also  propose  conditions  assuring  convergence  of  this  iterative  mapping 
to  a  unique  fixed  point.  We  further  compute  the  asymptotic  multiuser  effieieney  (AME)  [29]  of  this  overall 
system,  under  some  mild  assumptions  on  the  channel  decoders.  It  should  be  noted  that  the  analysis  in  this 
paper  is  based  on  large  sample  and  large  system  analysis. 

The  remainder  of  this  paper  is  organized  as  follows.  Section  II  introduces  the  signal  model  and  the  chan¬ 
nel  decoder  used  in  our  analysis.  The  performance  analyses  of  ML  channel  estimation  and  PIC  multiuser 
detection  are  given  in  Seetion  III  and  Seetion  IV,  respectively.  Based  on  these  results,  the  corresponding 
iterative  mapping  is  described  and  analyzed  in  Section  V.  Numerical  results  and  conclusions  are  given  in 
Section  VI  and  Section  VII,  respectively.  The  notations  used  in  this  paper  are  explained  as  follows. 

•  Throughout  this  paper,  if  no  special  note  is  given,  we  denote  vectors  with  small  letters  in  bold  fonts, 
matrices  with  capital  letters  in  bold  fonts  and  sealars  with  non-bold  fonts. 

•  For  any  variable  X,  we  denote  the  corresponding  estimate  from  the  decision  feedback  by  X  and  the 
corresponding  error  X  —  Xhy  5X. 
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•  Superscript  T  denotes  transposition  and  superscript  H  denotes  conjugate  transposition. 

•  I  denotes  the  identity  matrix. 

•  [x]  denotes  the  smallest  integer  larger  than  or  equal  to  x. 

•  mod(i,  j)  denotes  the  modulo  of  i  with  respect  to  j,  with  the  convention  of  mod(i,  i)  =  i. 

•  For  a  matrix  A^xn,  ||  A||p,  =  YTj=i  is  the  Frobenius  norm  of  A. 

II.  Signal  Model 

A.  Signal  Model 

We  consider  a  synchronous  uplink  long  code  (aperiodic)  DS-CDMA  system,  with  identical  channel  cod¬ 
ing,  binary  phase-shift  keying  (BPSK)  modulation,  K  active  users,  spreading  gain  A,  system  load  (3  =  ^, 
and  identical  transmission  rates  for  all  users.  The  transmitted  symbols  experience  multipath  fading.  We 
adopt  a  block  fading  model  and  denote  by  M  the  coherence  time,  measured  in  the  number  of  symbol  peri¬ 
ods,  over  which  the  channel  is  stationary.  Within  a  coherence  period,  the  chip  matched  filter  output  of  the 
receiver  at  symbol  period  t  can  be  collected  into  an  A-vector  given  by 

K  L 

v{t)  =  '^hk{t)'Y^akiSki{t) +  n{t),  t  =  l,2,...M,  (1) 

k=l  1=1 

where  L  denotes  the  number  of  resolvable  paths  per  user,  bk{t)  denotes  the  channel  coded  binary  symbols,  Oki 
denotes  the  channel  gain  of  the  /-th  path  of  user  k,  Ski{t)  denotes  the  binary  spreading  code  with  ||sfc/(f)  ||  =  1 
received  from  user  k  along  path  I  at  time  t  and  n(t)  is  an  A-vector  of  independent  and  identical  distributed 
(i.i.d.)  circularly  symmetric  complex  Gaussian  (CSCG)  ^noise  variables  with  (normalized)  variance  cr^.  It 
should  be  noted  that  although  the  assumption  of  synchronicity  is  valid  in  time  division  duplexing  (TDD) 
systems,  it  does  not  hold  for  many  frequency  division  duplexing  (FDD)  systems.  However,  as  it  will  be 
shown,  the  results  from  the  analysis  of  synchronous  systems  are  also  reasonably  valid,  though  not  exactly 
the  same,  in  the  case  of  asynchronous  systems. 

For  the  system  model,  we  have  the  following  assumptions. 

Assumption  II. 1:  The  channel  gains  {au}  are  independently  CSCG  distributed  with  zero  means  and  vari¬ 
ances  We  consider  only  the  case  of  large  L,  which  implies  that  'Yld=i  l^fczP  ~  1,  /c  =  1, ...,  K;  thus  all 
users  achieve  the  same  performance  with  maximal  ratio  combining  (MRC). 

complex  random  variable  is  CSCG  distributed  if  its  real  and  imaginary  parts  are  mutually  independent  Gaussian  random  variables  with 
zero  mean  and  identical  variance. 
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Assumption  II. 2:  We  ignore  intersymbol  interference  (ISI)  and  assume  that  the  spreading  codes  received 
along  different  paths  of  a  given  user  are  mutually  independent  {independent  model). 

Assumption  II. 3:  Based  on  Assumption  11.2,  the  crosscorrelations  Pkimn{l)  —  (note  that 

Pkiki{t)  =  1)  satisfy 

.  E  {pkimn{t)]  =  0,  if  {k,  1)  ^  (m,  n); 

•  E  [ph^S^)}  =  jf,  if  {k,  1)  ^  (m,  n); 

.  E  {pklmn{t)ppqrs{t)]  =  0,  if  {k,  /,  771,  u)  ^  (p,  g,  T,  s). 

The  above  assumptions  simplify  the  performance  analysis  substantially.  Moreover,  these  assumptions  are 
reasonable  for  practical  systems  due  to  the  following  reasons: 

•  Assumption  II.  1  is  based  on  the  fact  that  more  propagation  paths  are  resolvable  in  CDMA  systems 
than  narrow  band  systems,  particularly  in  environments  with  abundant  scattering  (e.g.,  indoor  environ¬ 
ment).  With  this  assumption,  we  ignore  the  impact  of  the  fluctuation  of  received  power  incurred  by  the 
multipath  fading,  and  consider  only  the  impairment  caused  by  the  channel  estimation  error. 

•  Assumption  II.2  is  unrealistic  since  these  sequences  are  shifted  versions  of  each  other  {shifted  model). 
However,  the  accuracy  of  the  results  dependent  upon  this  assumption  is  validated  with  numerical  results 
in  Section  VI  and  asymptotic  analysis  given  in  Appendix  1 . 

B.  Receiver  Structure 

The  structure  of  receiver  is  shown  in  Figure  1.  The  channel  coefficients  are  estimated  in  the  channel 
estimator,  which  operates  in  a  ‘semi-blind’  way.  Training  symbols  are  available  to  obtain  an  initial  estimate 
in  the  first  iteration.  In  the  further  iterations  the  information  symbol  decisions  from  channel  decoders  are 
assumed  to  be  correct.  Then,  both  the  training  symbols  and  fed  back  decisions  are  considered  as  training 
symbols  and  used  for  ML  channel  estimation.  A  multiuser  detector  is  used  to  mitigate  the  MAI  and  its 
outputs  are  de-interleaved  and  decoded  in  the  channel  decoder.  In  the  multiuser  detector,  we  use  the  LMMSE 
algorithm  in  the  first  iteration  and  the  PIC  algorithm  with  the  aid  of  hard  decision  feedback  in  the  succeeding 
iterations.  We  follow  the  standard  procedure  in  turbo  multiuser  detection  [1]  [13]  [22]  [30]  to  reconstruct  the 
channel  symbols  from  the  channel  decoder  output.  Then  these  channel  symbol  estimates  are  interleaved  and 
fed  back  to  the  multiuser  detector  and  channel  estimator  to  enhance  the  performance  iteratively. 

We  denote  by  bk{t)  the  estimated  binary  channel  symbol  of  user  k  at  symbol  period  t  that  is  fed  back  from 
the  channel  decoder.  For  simplicity,  we  use  hard  decision  feedback  and  denote  the  feedback  symbol  error 
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rate  by  P^.  The  decision  feedback  error  is  denoted  by  5hk{t)  =  bk{t)  —  bk{t).  Supposing  that  both  bk{t)  and 
5bk{t)  are  symmetrically  distributed,  it  is  easy  to  check  that 
.  E{6bk{t)}  =  0-, 

•  E  {bk{t)6bk{t)}  =  2Pe-, 

.  £«(«)}  =4P.. 

•  E  {6bk{m)6bi{n)}  =  0,  when  (/c,  m)  ^  (/,  n). 

It  should  be  noted  that,  in  practical  systems,  soft  decision  feedback  will  achieve  better  performance  than 
hard  decision  feedback.  However,  the  performance  of  channel  estimation  with  soft  decision  feedback  is 
determined  by  both  the  first  and  second  moments  of  the  decision  feedback  error  [17].  Thus  the  corresponding 
analysis  of  performance  evolution  is  more  complicated  than  the  case  of  hard  decision  feedback.  Therefore, 
we  adopt  hard  decision  feedback  in  order  to  simplify  the  system  performance  analysis. 

For  the  decision  feedback  from  channel  decoders,  we  have  the  following  reasonable  assumption,  which 
simplifies  the  analysis  and  is  also  used  in  [1]. 

Assumption  11.4:  The  codeword  length  is  assumed  to  be  large  enough  so  that  the  transmitted  symbols  are 
coded  over  many  coherence  periods.  The  decision  feedbacks  |6fc(t)  |  are  mutually  independent  for  different 
k  or  f. 


III.  Performance  Analysis  of  Channel  Estimation 

In  this  section,  we  discuss  the  performance  of  channel  estimation.  First,  we  explain  the  training  symbol 
based  ML  channel  estimation  algorithm  that  is  used  in  the  first  iteration.  Then,  we  consider  the  estimation 
of  the  channel  coefficients  with  only  hard  decision  feedback  from  the  channel  decoders.  Finally,  we  extend 
the  performance  results  to  channel  estimation  with  both  training  symbols  and  decision  feedback,  the  latter 
of  which  is  used  in  the  further  iterations. 

In  applying  the  turbo  principle,  to  avoid  the  reuse  of  information,  only  observations  are  used  in 

the  channel  estimation  for  multiuser  detection  in  symbol  period  i.  Thus  the  corresponding  channel  estimation 
error  is  independent  of  r(i).  However,  for  simplicity  of  discussion,  we  still  assume  that  all  M  received 
signals  are  used  for  the  channel  estimation  while  retaining  this  independence  assumption.  For  large  M,  this 
results  in  only  a  small  error  in  the  analysis. 

In  the  following  discussion  of  channel  estimation  and  PIC,  we  regard  the  channel  gains  {oki}  and  the 
spreading  codes  {s^i}  as  realizations  of  random  variables.  Only  the  transmitted  symbols,  decision  feedback 
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errors  and  noise  are  considered  as  random  variables.  Throughout  this  paper,  all  expectations,  denoted  as 
E{-},  are  over  the  distrihutions  of  these  three  variables.  Thus  our  results  are  conditioned  on  the  realizations 
of  {ttki}  and  {ski}.  However,  by  the  strong  law  of  large  numbers,  we  will  see  that  we  can  obtain  identical 
results  for  almost  every  realization  of  {aki}  and  {ski}  in  the  large  system  limit  {K,  N  oo). 

A.  Training  Symbol  Based  ML  Channel  Estimation 

First  we  assume  that  there  are  M  training  symbols,  channel  symbols  known  to  the  receiver,  within  a  single 
coherence  period.  For  simplicity  in  deriving  the  channel  estimate,  we  stack  the  chip  matched  filter  output  of 
the  signal  corresponding  to  these  training  symbols,  rewriting  (1)  as 

r  =  Sa  +  n,  (2) 


where 


n 

a 

S 


H 

AfMxl  ’ 
H 

NMxl 


—  (oii,  ai2, okl) 


KLxl  > 


=  (S(l)B(l)f  ,...,(S(M)B(M))^ 


NMxKL 


( 


B(m)  = 


bi{m)lLxL  0 

0  b2(m)lLxL 


0 

0 


0 


\  0  0  ■■■  bKim)lLxL  I 

S(m)  =  (sii(m),Si2(m),...,Sifi(m))^^^^,  m  =  l,...,M. 

Applying  the  ML  criterion  and  the  normality  of  the  noise,  we  can  obtain  the  ML  channel  estimate,  which 
is  given  by 

a  =  argmaxP(r|a) 

a 

=  argmin  ||r  —  Sa|| 

a 

=  (S^S)”iS^r 

=  R'V,  (3) 
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where  R  =  S^S  and  y  =  S^r. 

It  follows  direetly  that  the  ehannel  estimation  error  is 


—  cl  -  cl 

=  — R  ^S^n, 


from  which  it  is  obvious  that  this  error  has  zero  mean  and  covariance  Sa  =  -E  {(5a(5a^}  =  R  \ 

For  a  finite  M,  we  can  compute  trace  {R“^ }  in  the  large  system  limit  (i.e.  when  K,  N  ^  oo  while  keeping 
the  system  load,  ^  =  (3,  constant).  For  a  system  with  system  load  (3,  it  is  well  known  that  as  K  ^  oo. 


trace{R-i}  to  the  multiuser  efficiency  of  a  decorrelator,  namely  I  —  (3  [29].  is  equivalent  to  the 

livalent  sysh 

trace  {Sa} 


covariance  matrix  of  a  system  with  equivalent  system  load  (3'  =  =  -^(3.  Thus  as  K,N  oo,wc  have 


M 


M  -  L(3' 


Therefore,  for  sufficiently  large  K  and  N ,  the  variance  of  channel  estimation  error  is  given  by 

.2 


A„  = 


O’; 


M-Lp' 


(4) 


which  can  be  approximated  by  ^  when  M  is  sufficiently  large. 

It  should  be  noted  that,  in  asynchronous  systems,  we  can  remove  part  of  the  chips  in  the  first  and  the 
last  symbol  periods  to  obtain  a  similar  matrix  SNM-dm,axxKL,  where  d^ax  denotes  the  largest  time  offsets 
of  different  users,  measured  in  chips.  Since  the  training  symbols  have  been  incorporated  into  the  spreading 
codes,  we  can  consider  the  columns  of  S  as  random  {NM  —  dmax)—  vectors,  regardless  of  the  time  offsets 
of  different  users.  Therefore,  the  variance  of  channel  estimation  error  in  asynchronous  systems  is  similar  to 
that  of  synchronous  systems  when  M  is  sufficiently  large. 


B.  Channel  Estimation  with  Decision  Feedback 

1)  Algorithm:  When  decision  feedback  is  used  in  place  of  training  symbols  to  derive  the  ‘ML’  channel 
estimates^,  a  process  that  assumes  that  the  decision  feedback  is  free  of  error,  the  channel  estimation  error 
is  caused  by  both  the  thermal  noise  and  the  decision  feedback  error.  On  applying  (3),  the  channel  estimate 

^By  ‘ML’  estimates,  we  mean  using  the  expression  obtained  from  the  training  symbol  based  estimation,  but  with  symbols  obtained  from 
decision  feedback.  It  is  not  an  exact  ML  estimate  since  the  distribution  of  the  decision  feedback  error  is  not  considered. 
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with  decision  feedback  is  given  by 


a  =  R“V 

=  R-^S^(Sa  +  n) 

=  a  +  R~^S^((5Sa  +  n), 

where  =  S  -  S,  R  =  S^S,  y  =  S^r  and  S  is  the  version  of  S  in  (3)  obtained  from  the  decision  feedback, 
which  is  given  by 

S  =  f(s(l)B(l))^,...,(s(M)B(M))^') 

^  ^  '  / NMxKL 

bi{m)lLxL  0  ■■■  0 

0  b2{m)lLxL  ■■■  0 

0  0  ■■■  bK{m)lLxL 

Hence,  the  channel  estimation  error  can  be  decomposed  into  two  parts 

(5a  =  -R-iS^((5Sa  +  n) 

=  <5aj  -f  (5ajj,  (5) 

where  (5a/  =  — R“^S^(5Sa  and  (5a„  =  — R^^S^n  denote  the  channel  estimation  error  due  to  the  decision 
feedback  error  and  the  thermal  noise,  respectively.  It  is  reasonable  to  assume  that  (5a/  and  (5a„  are  mutually 
independent.  (Recall  our  assumption  concerning  the  use  of  only  measurements  t  ^  i  in  estimating  gains  at 
time  i.) 

It  is  difficult  to  tackle  the  calculation  of  (5a  due  to  the  matrix  inversion  R^^.  However,  we  can  approximate 
R“^  by  when  Pe  is  sufficiently  small.  This  approximation  is  justified  by  the  following  lemma. 

Lemma  II LI:  When  fixing  K  and  N,  we  have 

M'R~^  Iklxkl, 

almost  surely  ^  as  M  — >  cx)  and  Pg  — ^  0. 

Proof:  According  to  the  definition  of  R,  we  have 

R-i  =  R  1  +  R  lA, 

^Here,  a  matrix  is  considered  as  a  point  in  the  probability  space  and  the  metric  is  induced  by  a  matrix  norm. 
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where  A  =  (I  —  (5RR 

E{\\A\\f}<E 

whieh  tends  to  0  as  Pg  — ^  0.  Thus,  we  have 


I.  Aecording  to  the  error  analysis  of  matrix  inversion  in  [1 1],  we  have"^ 

II^RR-'IIf 


1  -  ||5RR^1| 


=  0(^e), 


R~^  -  R 


<  R-i  ,P{||A||^}^0, 


as  Pg  — >  0.  Therefore,  R  ^  converges  to  R  ^  almost  surely  as  Pg  — >  0. 

Applying  the  strong  law  of  large  numbers  and  the  fact  that  the  diagonal  elements  in 


M 


(B(m)S(m))'^S(m)B(m) 


m=l 


are  M  and  the  off-diagonal  elements  in  (B(m)S(m))^S(m)B(m)  are  independent  for  different  values  of  m 
and  have  zero  mean,  we  obtain  that,  while  keeping  K  and  N  fixed,  ^  ^  ^klxkl  almost  surely,  as  M  — cx). 
Since  the  elements  of  R“^  are  continuous  functions  of  those  in  R  in  a  neighborhood  of  R  =  MIklxkl,  we 
also  have  MR~^  — ^  ^klxkl  as  M  — >  cx).  This  completes  the  proof.  ■ 

Therefore,  we  can  further  approximate  R~^  by  for  large  M  and  small  Pg.  For  simplicity,  our 

further  discussion  of  (5a/  will  be  based  on  this  approximation,  which  will  be  validated  by  numerical  results. 
Consequently,  in  the  following  discussions,  we  use  the  approximations 


(5a/  =  - 


and 


2)  Covariance  matrix  of  channel  estimation  error:  We  denote  the  covariance  matrices  of  (5a,  5a/  and 
5a„  by  Sa,  Sf  and  Sn,  respectively,  which  satisfy  Sa  =  Sf  +  Sn.  We  first  consider  the  channel  estimation 
error  incurred  by  decision  feedback  errors.  The  following  lemma  shows  that  the  channel  estimation  error 
5a/  is  asymptotically  biased.  The  proof  is  given  in  Appendix  II. 

Lemma  III. 2:  When  keeping  K  and  N  fixed,  we  have 


P{5a/}  ^  2Pga,  (6) 

almost  surely,  as  M  — >  cx. 

*x  =  0{Pe)  means  <  oo  as  Pe  ^  0. 
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It  should  be  noted  that  this  bias  cannot  be  removed  a  priori  in  the  estimator  since  it  is  dependent  on  the 
channel  gain,  a.  However,  this  bias  vanishes  as  Pg  — ^  0. 

An  asymptotic  expression  for  the  elements  in  Sf  is  given  in  the  following  proposition,  whose  proof  is 
given  in  Appendix  III,  where  we  also  explain  that  the  conclusion  also  applies  to  asynchronous  case  when  Pg 
is  sufficiently  small. 

Proposition  III.3:  For  all  i  and  j,  when  fixing  K  and  N,  we  have  that  (recall  that  Oki  is  the  channel  gain 
of  use  k  and  path  /) 


4Pe  Ori],  modiip)?  +  J:  Sfil, 

,  mod(fc,L)  P)  >  if*—  7) 

A/x(S,)y  - 

< 

4Pe  (l  +  ]v)  mod{j,Ly 

if*  ^7  and  [f]  =  , 

(V) 

^  4Pg  (l  +  mod{j,Ly 

if  Til  ^  Til 

almost  surely,  as  M  ^  oo. 

For  5a„,  which  is  caused  by  thermal  noise,  the  corresponding  analysis  is  identical  to  that  of  training 
symbol  based  estimation.  Then,  we  have 

MSn  =  Mcov 

= 

— >  O'n^KLxKL,  (8) 


almost  surely,  as  M  ^  oo.  Then  the  covariance  matrix  of  channel  estimation  error  Sa  =  P  {(5a(5a^}  = 
Sf  +  Sn  can  be  obtained  from  (7)  and  (8). 

3)  Variance  of  channel  estimation  error:  The  variance  of  channel  estimation  error  can  be  obtained  as  a 
corollary  of  the  previous  subsection. 

Corollary  IIL4:  On  defining  =  ^trace  {Sa},  we  have 


4Pg(l  +  /3L)  ,  2 

- F - ^ 


(9) 


almost  surely,  as  K,  N,  M  ^  oo. 

Thus,  when  K,  N,  M  are  sufficiently  large,  we  have  the  following  approximation 


A, 


4Pg(l  +  /3L) 

LM  M' 


(10) 


It  should  be  noted  that  the  channel  estimation  error  cannot  be  removed  by  increasing  M  although  the 
variance  vanishes  as  M  — >  cx),  since  the  estimate  is  biased  and  the  bias  cannot  be  removed  a  priori. 
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C.  Estimation  with  Both  Training  Symbols  and  Decision  Feedback 

We  denote  the  number  of  training  symbols  by  Mt  and  the  corresponding  percentage  by  a  =  ^ .  When  the 
training  symbols  and  decision  feedback  are  combined  for  channel  estimation,  the  performance  is  determined 
by  (10),  with  Pe  replaced  by  (1  —  a)Pe-  Decision  feedback  should  only  be  used  along  with  the  training 
symbols  if  the  resulting  variance  is  smaller  than  that  obtained  when  only  the  training  symbols  are  used. 
Then  it  is  easy  to  check  that,  when  M  and  Mt  are  sufficiently  large,  Pemax^  the  maximum  P^.  assuring 
performance  improvement  when  decision  feedback  is  used,  is  determined  by 

A{1  -  a)P^{l  +  I3L)  ^  ^ 

LM  M  -  Mt 

which  results  in 


p 

4a(l  +  /3L)’ 


(11) 


from  which  we  observe  that  Pemax  decreases  with  a  and  while  increasing  with  cr^  and  L. 


IV.  PIC  AND  Channel  Decoder 

A.  Performance  Analysis  of  PIC 

For  convenience  of  analysis,  the  performance  of  PIC  is  analyzed  based  on  matched  filter  (MF)  outputs. 
We  drop  the  index  of  the  symbol  period  for  notational  simplicity  throughout  this  section.  For  a  given  symbol 
period,  the  MF  outputs,  which  form  sufficient  statistics  for  multiuser  detection,  are  given  by 

y  =  S'^r. 


In  PIC  based  multiuser  detection,  the  MAI  reconstructed  from  the  channel  estimates  and  the  decoder  output 
is  subtracted  directly  from  the  MF  output  of  the  desired  user.  Without  loss  of  generality,  we  take  the  /-th 
path  of  user  1  as  an  example;  then  the  MF  output  after  PIC,  which  is  contaminated  by  residual  MAI  and 
thermal  noise  nu  =  s^^n,  is  given  by 


where 


Vll  —  ttilbi  +  dlmPlllmbl  +  T 
mp 


Ih 


K  L 


hi  —  EE  Pllkm  (  (^kmbk  )  P 

k=2  m=l 


(12) 
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which  is  the  sum  of  the  residual  interference  and  the  thermal  noise.  It  is  obvious  that  E{Iii}  =  0.  And  the 
corresponding  variance  is  given  by 


cr]^E{\Iu\^}  = 

k=2  m=l 
K  L 

=  {\Sakmf}  +  4Pe  \akmf  +  ^PeE  {\Sakmf}  +  2E  {Sa^n^aU  E  {bM 

k=2  m=l 

-2E  [\5akm\^]  E  {hkbhk}  -  8PeE  {akmSal^}  |  +  cr^ 

— i>  /3LAa  +  4/3(1  —  Pe)Pe  +  O'n)  (13) 


as  K,N  cx),  where  we  have  applied  the  faet  that  E  {\Sakmf}  =  +  4P|  \akmf,  E  {akmbal^}  = 

2Pe  \akmf,  E  {bkSbk}  =  2Pe.  It  is  easy  to  check  that  a]  is  identical  for  asynchronous  systems  since  different 
time  offsets  do  not  affect  the  interference  power. 

It  is  diffleult  to  apply  the  central  limit  theorem  to  show  the  asymptotic  normality  of  the  PIC  output  sinee 
the  variables  {5akm}  are  mutually  correlated  across  different  users  and  paths.  However,  numerical  results  in 
Section  VI  will  show  that  the  output  distribution  of  PIC  can  be  well  approximated  by  a  Gaussian  distribution. 
Thus,  in  the  subsequent  sections,  we  assume  that  the  output  of  PIC  is  Gaussian  distributed. 

According  to  the  properties  of  the  crosscorrelation  given  in  Section  II.A,  pmm  — ^  0  almost  surely,  as 
A  — CX).  Thus,  for  large  spreading  gain,  the  interference  across  different  paths  of  the  same  user  can  be 
ignored.  With  the  normality  assumption  of  the  residual  MAI,  it  is  easy  to  show  that  the  variables  ^ 

are  mutually  independent  as  N  ^  oo,  whieh  means  that  ehannel  coded  symbol  bi  is  transmitted  through  L 
independent  channels.  This  assumption  simplifies  the  analysis  although  it  does  not  hold  exactly  when  N  is 
finite.  Thus,  we  use  MRC  to  colleet  these  L  replicas,  resulting  in  the  output 

L  L 

^  dhaubi  +  ^  dhlii.  (14) 

i=i  i=i 


Applying  Lemma  III.2,  we  obtain  that,  as  M,  L  — >  cx, 

L  L 

'^^dhaiibi  =  {\aiif  —  Sahaii) 


1=1 


1=1 


1  -  y/  B  {ia;,} 


OlZ 


1=1 


i-2f.E 


On 


1=1 


=  1-2P, 
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Moreover,  we  can  obtain  that,  as  M,  L  ^  oo 


E 


1=1 


1=1 


^  1=1  1=1 
^  {l-AP,  +  AP^  +  LAa)aj 

=  ((l-2Per  +  LA>2. 


Therefore,  when  M  and  L  are  sufficiently  large,  (14)  can  be  approximated  by 


a 


2 

I 


Zi  ^  {1  -  2Pe)bi  +  ni,  (15) 

where  rii  is  a  CSCG  random  variable  with  variance  of  ((1  —  2Pe)^  +  LAajaj.  An  interesting  observation 
is  that  the  channel  estimation  error  not  only  increases  the  interference  but  also  decreases  the  valid  received 
power  of  the  desired  user. 


B.  Performance  of  Channel  Decoder 

At  the  channel  decoder,  Pe  is  a  function  of  the  input  signal-to-interference-plus-noise  ratio  (SINR)  at  the 
input  to  the  channel  decoder  given  by 

^  (sink)  ’ 

where  the  function  g  can  be  estimated  using  Monte  Carlo  simulations.  For  most  practical  channel  codes,  the 
following  assumption  is  reasonable: 

Assumption  IV.l:  Within  a  closed  interval  =  [0,  function  g  satisfies 

•  g{x)  monotonically  increases  with  x,  and  5f(0)  =  0; 

•  g{x)  is  continuously  differentiable  and  (^'(O)  =  0. 

V.  Analysis  of  System  Performance 

In  this  section,  we  analyze  the  overall  iterative  system  shown  in  Figure  1.  We  consider  only  the  case 
of  small  Pe,  moderate  and  moderate  M  and  note  that  the  analytic  results  become  more  precise  as  Pg 
and  (7^  decrease  and  M  increases.  This  configuration  is  reasonable  for  the  decision  feedback  basedsystems 
since  if  M  is  large,  training  symbol  based  channel  estimation  can  be  adopted  with  marginal  loss  of  spectral 
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efficiency;  if  M  is  small,  it  is  difficult  to  carry  out  coherent  detection;  and  if  Pg  is  large,  the  iteration  diverges. 
Although  the  performance  analysis  of  the  channel  estimation  in  Section  III  is  based  on  large  M,  numerical 
results  in  Section  VI  indicate  that  expression  (10)  is  still  valid  for  moderate  M.  We  adopt  the  expressions 
(10)  and  (13)  in  large  system  limits  {K,  N  oo). 

A.  Iterative  Mapping 

In  this  section,  we  consider  the  d-th  iteration  and  couple  the  results  from  Section  III  and  Section  IV 
to  analyze  the  overall  system  performance.  We  can  regard  the  decoding  process  as  an  iterative  mapping 
:  M  — M  in  terms  of  the  error  probability  of  the  decoder  output  after  the  d-th  iteration,  ,  which  is 
given  by  (recall  that  g  is  defined  as  the  function  characterizing  the  output  error  probability  in  terms  in  input 
SINK  in  (16)) 

pW  =  d(Pi'^-i)) 

^  (Po  +  ,  (17) 

where  we  ignore  terms  of  a  smaller  order  than  Pg  and  P  since  we  assume  small  Pg  and  large  (or  moderate) 
M.  Based  on  (10),  (13)  and  (15),  the  coefficients  Dq  and  Di  are  given  by 

Po  =  o-n  (l  +  fy  + 

Pi  -  4  I ^ - 

B.  Condition  for  Convergence 
A  reasonably  good  initialization,  which  results  in  sufficiently  small  channel  estimation  error  and  MAI 

in  the  first  iteration,  is  necessary  to  guarantee  the  convergence  of  the  iterative  mapping  described  in  (17). 
In  the  initial  stage,  only  training  symbols  are  used  for  the  channel  estimation  since  no  decision  feedback 
is  available  then.  Any  non-iterative  multiuser  detection  technique  can  be  applied  to  the  initializing  stage. 
For  practical  applications,  we  can  use  the  LMMSE  detector,  whose  performance  using  imperfect  channel 
estimation  can  be  obtained  using  the  replica  method  [18]. 

For  convergence,  the  variance  of  input  interference  and  noise  of  the  initializing  stage,  denoted  by  aj{0) 
and  obtained  from  the  SINK  of  the  LMMSE  detector,  must  satisfy  the  following  conditions: 

•  cr|(0)  is  located  within  the  interval  P  defined  in  Section  IV.B,  namely 

a^O)  <  (18) 
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This  condition  assures  a  reasonably  good  initial  performance  of  the  iterations. 
•  The  variance  of  interference  and  noise  decreases  with  iteration  time,  namely 

aj{0)  -  Do 


Di 


This  condition  assures  that  the  iterations  do  not  diverge. 


(19) 


C.  Condition  Assuring  the  Uniqueness  of  the  Fixed  Point 

If  there  exists  more  than  one  fixed  point,  the  iteration  may  beeome  stuck  at  a  suboptimal  fixed  point  and 
not  converge  to  the  optimal  one.  The  following  proposition  provides  a  sufficient  condition  for  the  uniqueness 
of  the  fixed  point  and  the  eorresponding  convergence  rate. 

Proposition  V.l:  (1)  If  there  exists  a  7  <  1,  such  that 


max^-gQ  {g'{x)) 

then  there  exists  only  one  fixed  point  Xf  for  the  iterative  mapping  Xk+i  =  h  (xk),  and  for  every  initial  point 

/c 

Xq  G  f),  the  mapping  eonverges  to  Xf  with  an  exponential  rate,  namely  \\xk  —  Xf\\  <  ||a;o  —  Xf\\. 

(2)  If  there  exists  an  xi  G  O  such  that  <  Di  <  then  there  exists  a  Dq  such  that  there  is  more 
than  one  fixed  point  for  h. 

Proof:  (1)  The  condition  Di  <  max  J^{g'{x))  implies  that  h!{x)  =  g'{Do  +  Dix)  <  7  <  1.  Then  hf) 
is  a  contraction  mapping,  and  the  conclusions  follow  due  to  Banach’s  fixed  point  theorem  [14]. 

(2)  Letting  Xf  =  g{xi)  and  setting  Do  =  Xi  —  DiXf,  we  can  show  that  Dq  >  0  due  to  the  assumption  that 
Di  <  It  is  easy  to  check  that  Xf  is  a  fixed  point  and  g'{Do  +  DiXf)  =  Digfxi)  >  1.  Hence, 

there  exists  an  e  >  0  such  that  for  all  x  G  ixf,Xf  +  e),  g{Do  +  Dix)  >  x.  However,  g{Do  +  D1X2)  <  X2 
for  X2  =  g  (o'KO))  due  to  condition  (19).  If  X2  <  Xf,  there  exists  at  least  one  fixed  point  within  (0,  X2)  since 
g{Do)  >  0;  if  X2  >  Xf,  there  exists  at  least  one  fixed  point  different  from  Xf  within  {xf,  X2).  ■ 

It  should  be  noted  that  condition  (20)  is  sufficient  but  not  necessary  for  the  uniqueness  of  the  fixed  point. 
This  condition  is  more  stringent  than  the  condition  of  convergence  in  (19)  since  it  assures  both  the  uniqueness 
of  the  fixed  point  and  the  exponential  eonvergence  rate.  The  second  part  shows  that  a  moderate  Di  may  cause 
multiple  fixed  points.  A  useful  conclusion  drawn  from  (20)  is  that  this  iterative  procedure  does  not  work 
well  for  those  channel  codes,  such  as  powerful  turbo  eodes  or  LDPC  codes,  that  have  a  steep  performance 
curve  (bit  error  rate  versus  SINK)  which  implies  a  large  value  of  maxa-g^  (gfx)).  This  will  be  demonstrated 
in  numerical  simulations  in  Section  VI. 
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D.  Asymptotic  Multiuser  Efficiency 

As  is  described  in  [29],  the  asymptotic  multiuser  efficiency  measures  the  slope  at  which  the  bit-error-rate 
goes  to  zero  in  logarithmic  scale,  giving  intuition  into  the  performance  loss  from  multiuser  interference. 

Suppose  that  there  is  only  one  fixed  point  for  the  iterative  mapping  h,  and  let  Pe(o'n)  be  this  fixed  point 
when  the  noise  power  is  Similarly,  let  i9o(a^)  and  Di{a'^)  be  the  corresponding  values  of  Dq  and  Di  in 
(17).  It  is  obvious  that  Pe(0)  =  0  and  /i>o(0)  =  0. 

The  asymptotic  multiuser  efficiency  is  given  by 


AME  =  lim 


cr„ 


.1^0  Doiaf)  +  D,ial)P,ial) 
1 


dPoPl) 


dal 


tI=0 


+ 


d(DPal)P,(al)) 


dal 


’-1=0 


If  H{Pe,  af)  =  g  (E>o(<)  +  Di{af)Pe)-Pe,  then  Peiaf)  is  the  unique  solution  of  i7(Pe,  <)  =  0-  Applying 
the  assumptions  that  (^'(O)  =  0  and  Pe(0)  =  0,  we  have 

d{D,{al)PMl)) 


daf 


>■2=0  da?, 


>-2=0 


=  -fli(O) 

=  -fl.(o)- 
=  0. 


da^ 


0-2=0 


dH(P,,al) 

dPe 

’’1=0 

d{Do{af)+Dl{al)Pe) 


dal 


’1=0 


g\^) 


D^{f))g'{f)  -  1 


Thus 


AME  = 


dPoial)  I  , 

dal  \<=0 


1  ,  M' 

M 


(21) 


Erom  (21),  we  can  see  that  the  loss  of  AME  is  due  to  the  channel  estimation  error  incurred  by  the  thermal 
noise.  The  impact  of  the  decision  feedback  error  vanishes  as  af  0,  while  that  of  the  channel  estimation 
error  remains. 


E.  Computational  Aspect 

The  main  computational  cost  of  the  iterative  channel  estimation  and  multiuser  detection  includes: 
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•  Solving  the  linear  equation  Ra  =  y  for  ML  ehannel  estimation. 

•  Reconstrueting  the  channel  symbols  and  cancelling  the  interference. 

•  Channel  decoding. 

Since  the  channel  symbol  reconstruction  is  similar  to  the  encoding  procedure  and  the  interference  cance¬ 
lation  requires  only  subtractions,  this  is  not  a  bottleneck  of  the  whole  procedure  and  the  corresponding 
computational  cost  is  of  complexity  0{K).  Real-time  channel  decoding  can  also  be  accomplished  in  a  way 
similar  to  Turbo  codes.  Therefore,  the  main  bottleneck  is  solving  the  linear  equation  for  channel  estimation. 

Direct  Gaussian  Eliminatation,  which  is  of  complexity  0{K^),  can  be  applied  to  solve  the  equation  Ra  = 
y  when  K  is  small.  When  K  is  large,  iterative  techniques  of  solving  linear  equations,  such  as  the  Jacobi 
method  and  the  Gauss-Seidel  method,  can  be  applied.  For  assuring  the  convergence,  we  cite  the  following 
lemma  from  [10]: 

Lemma  V.2:  The  sufficient  and  necessary  condition  for  the  convergence  of  iterations  in  solving  the  linear 
equation  Ax  =  y  is  that 

•  A  and  2  diag(A)  —  A  are  both  positive  definite  in  the  Jacobi  method^; 

•  A  is  positive  definite  in  the  Gauss-Seidel  method. 

The  Gauss-Seidel  method  always  converges  when  (3  <  1  since  R  is  positive  definite  when  K  <  N.  For 
the  Jacobi  method,  it  is  easy  to  check  that  diag(R)  =  \kxk-  Since  the  largest  eigenvalue  of  R  converges  to 
(l  -|-  a/^)^  [3]  almost  surely  as  K,N  oo,  the  eigenvalues  in  2  diag(R)  —  R  are  less  than  2  —  (l  -|-  a/3)^ 
almost  surely  in  the  large  system  limit.  Therefore,  <  1  is  a  sufficient  condition  for  the  almost  sure 
convergence  of  Jacobi  iteration  in  the  large  system  limit.  Then,  when  K  and  N  are  sufficiently  large  and 
K  <  N,  we  can  use  either  Gauss-Seidel  or  Jacobi  iterations  to  estimate  the  channel  coefficients  efficiently. 

VI.  Numerical  Results 

A.  Channel  Estimation 

Figure  2  shows  the  average  variance  of  the  channel  estimates  versus  the  coherence  time  M  with  the 
configuration  of  j3  =  0.2,  L  =  5,  =  0,  Pg  =  0.1  and  the  signal-to-noise  ratio  (SNR)=  5dB^.  The 

asymptotic  results  obtained  from  (10)  and  the  simulation  results  for  finite  systems  {N  =  100)  with  spreading 
codes  for  the  shifted  model  are  represented  by  solid  and  dotted  curves,  respectively.  In  this  figure,  the 

®diag(X)  denotes  a  diagonal  matrix  constituted  by  the  diagonal  elements  in  matrix  X 

®Note  that  Pe  and  SNR  are  not  mutually  independent;  however,  we  set  these  two  parameters  arbitrarily  to  test  the  validity  of  asymptotic  results. 
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estimation  error  varianee  eaused  by  deeision  feedback  and  noise  are  denoted  by  A/  and  A„,  respectively. 
The  corresponding  asymptotic  results  are  obtained  from  the  first  and  the  second  terms  in  (10),  respectively. 
We  can  observe  that  the  asymptotic  results  match  the  simulation  results  well  even  when  M  is  small.  This 
figure  also  demonstrates  the  validity  of  results  based  on  the  independence  assumption  of  the  spreading  codes 
given  in  Section  II.A. 

B.  Normality  of  PIC  Output 

Figure  3  shows  the  channel  symbol  error  rate’  with  the  configuration  of  SNR  =  lOdB,  A  =  iV  =  30 
and  Pe  =  0.1  and  0.05.  The  solid  curves  represent  the  results  obtained  from  numerical  simulations  and  the 
dashed  curves  represent  the  results  with  the  assumption  that  the  output  of  PIC  is  CSCG  distributed.  The  gap 
between  the  numerical  results  and  CSCG  based  prediction  is  small,  thus  justifying  the  normality  assumption 
of  the  PIC  output. 

C.  User  Capacity 

We  define  the  user  capacity  to  be  the  maximum  system  load  Pmax  with  which  the  system  can  achieve 
the  information  bit  error  rate  of  10“^.  Two  types  of  channel  codes,  the  convolutional  code  (35,  23)8  a 
turbo  code  (with  two  constituent  codes  (37,  21)3),  with  bit  rate  R  =  \  and  codeword  length  1024  are  used 
in  this  paper  and  their  error  rates  for  both  information  bits  and  extrinsic  information  based  channel  symbols 
are  shown  in  Figure  4.  The  corresponding  /Umax’s  for  various  values  of  coherence  time  M,  denoted  by 
‘iterative’,  are  given  in  Figure  5  and  Figure  6  for  convolutional  codes  and  turbo  codes,  respectively,  with  the 
configuration  a  =  0.2,  SNR=  5dB  and  L  =  5.  The  (3max'^  of  the  non-iterative  LMMSE  detector,  denoted 
by  ‘LMMSE’,  are  given  for  comparison.  We  can  see  that  the  iterative  system  achieves  substantially  higher 
user  capacity  than  the  non-iterative  one.  The  performance  of  systems  with  ideal  initialization,  where  actual 
channel  parameters  are  provided  by  a  genie  in  the  initialization  stage,  denoted  by  ‘Perfect  initialization’, 
implies  that  a  good  initialization  can  improve  the  performance  considerably.  Thus,  blind  or  semi-blind  non¬ 
iterative  techniques,  which  make  use  of  information  symbols,  can  be  applied  to  obtain  a  better  initialization. 
Eor  comparison,  the  user  capacities  of  both  iterative  and  non-iterative  systems  with  perfect  channel  state 
information  are  also  given  in  both  figures.  An  interesting  observation  is  that  the  relative  performance  gain 

^This  channel  symbol  error  rate  is  equivalent  to  bit  error  rate  when  the  output  of  PIC  is  used  directly  for  the  detection  (without  channel 
decoding). 
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of  iterative  systems  over  the  non-iterative  ones  is  smaller  for  turbo  codes  than  for  convolutional  codes.  This 
is  due  to  the  steeper  waterfall  region  in  turbo  codes. 

VII.  Conclusions 

In  this  paper,  we  have  analyzed  the  performance  of  decision  feedback  based  iterative  channel  estimation 
and  multiuser  detection  in  multipath  DS-CDMA  channels.  The  decoding  process  has  been  described  as  an 
iterative  mapping  in  terms  of  the  variance  of  the  channel  decoder  output,  and  conditions  assuring  the  conver¬ 
gence  and  uniqueness  of  a  fixed  point  have  been  proposed.  Numerical  results  show  that  the  initialization  is 
important  to  the  iterations,  thus  necessitating  the  use  of  non-iterative  blind  or  semi-blind  channel  estimation 
algorithms  for  initialization  purposes.  Another  observation  of  interest  is  that  the  gain  of  the  iterative  process 
over  a  non-iterative  one  is  small  when  a  near-optimal  channel  coding  scheme  is  used. 

Appendix  I 

Validity  of  Independence  Model  for  Spreading  Codes 

In  (1),  for  different  values  of  I  and  m,  s^i  and  are  generated  by  the  same  binary  sequence  with  different 
offsets.  Our  purpose  is  to  show  that  if  K  and  N  are  large  enough,  we  can  regard  the  shifted  spreading  codes 
of  different  paths  of  a  given  user  as  independent  sequences.  The  properties  based  on  this  assumption,  which 
are  used  for  the  system  performance  analysis  in  this  paper,  include: 

•  The  properties  of  crosscorrelation  pkimn  in  Section  II.A. 

•  The  distribution  of  the  eigenvalues  of  the  matrix  SS^,  when  developing  the  expression  of  A„  for  finite 
M  and  large  K  in  Section  III.C.  Our  assumption  means  that  the  corresponding  distribution  of  the  shifted 
model  is  asymptotically  identical  to  that  of  the  independent  model. 

It  is  easy  to  check  the  first  item  using  the  symmetry  of  the  binary  distribution.  However,  the  validity  of 
the  second  one  is  non-trivial  and  is  of  considerable  importance  when  applying  the  theory  of  large  random 
matrices  to  multipath  fading  channels.  We  can  tackle  this  problem  by  showing  that  the  moments  of  the 
eigenvalues  in  both  models  are  the  same  via  the  following  lemma. 

Lemma  I.l:  Denote  a  generic  eigenvalue  of  SS^  by  A.  Then  the  m-th  moment  of  A  in  the  shifted  model 
is  given  by 

m 

E  =  'Y^{j3'f  ^  c(mi,  ...,mfc),  as  K  ^  oo, 

k=l  mi+...+m;c=m 
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which  is  the  same  expression  of  that  of  the  independent  model,  and  where  the  definition  of  c(mi,  ...,mk)  is 
given  in  [19]  and 

Proof:  Using  similar  arguments  to  those  in  [19],  we  have 


-B{trace{(SS’')”'}} 


K 


N 


Y.  Y.  E[v, 


^1  — 1  dm  }  7 


(22) 


where  Ujj  =  V^S^. 

For  any  v  7^  is,  Vir,jp  =  ^is,jq  when  \f~\  =  \f~\  and  jp  —  jq  equals  the  offset  difference  between  these 
two  shifted  sequences.  However,  the  probability  of  such  events  vanishes  as  K  ^  oo  sinee 


/  rn  \  2L  +  1 

P(|v -i^l  <  L)  <  (  j  - ' asK^oo. 

Thus,  as  K  ^  oo,  the  term  involving  l^j  ’s  of  different  users,  which  are  mutually  independent,  dominates 
the  summation  in  (22).  The  remaining  part  of  the  proof  is  the  same  as  in  [19].  ■ 

The  following  lemma  (Theorem  30.1  in  [5])  provides  a  sufficient  eondition  for  the  equality  of  two  proba¬ 
bility  measures  when  their  moments  are  identieal . 

Lemma  1.2:  Let  /i  be  a  probability  measure  on  the  real  line  having  finite  moments  ak  =  x^fi{dx)  of 
all  orders.  If  the  power  series  has  a  positive  radius  of  eonvergenee,  then  jj,  is  the  only  probability 

measure  with  the  moments  {om}m=i  2  • 

For  applying  Lemma  1.2,  we  need  the  following  lemma  which  provides  an  upper  bound  for  the  moments 
of  the  eigenvalues. 

Lemma  1.3:  For  any  eigenvalue  A  of  SS^,  there  exists  a  eonstant  C  >  max(l,/)')  such  that  for  m  = 

1,2,... 


E  {\^}  <  .  (23) 

Proof:  The  result  follows  by  induetion  on  m. 

It  is  easy  to  verify  that  (23)  holds  when  m  =  1,  2.  Suppose  E  {A*^}  <  for  n  =  1,  2, ...,  m.  Use 

the  following  recursive  formula  [19]  to  evaluate  E  whieh  is  given  by 

m+1 

E{\^+^]  =  ^I3'  . 

k=l  mi+...+m;c=m+l 
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Then  we  have 


m— 1 


E{X^+^}  =  p’ [l  +  mE{X}  +  E{X^}  +  J2  Y1  E  {X^^-^}  ...E  {X^^-^} 

k=2  mi+...+mfc=mH-l 
m—1  k 


<  (3'  {l  +  m(3'  +  ^  ^  JJ 

\  k=2  mi+...+mfc=mH-l  2=1 

(m—1 

l  +  mp'  +  C^rrT-^  +  ^  ^  C 

k=2  mi+...+mj.=m+l 

/ 

m—1 


mj— 3 


m 


m+l— /c^m— 1— /c 


<  C'”^+1  1  +  ^ 


m—2 


<  C'”^+1  1  +  ^ 


1  m  \ 

1 

U-J 

J 

/  m  —  1  ^ 

\ 

\  k 

1  / 

,  ^  - 1 

= 

k=o  y  k 


m 


m—l—k 


where  the  first  inequality  is  based  the  assumption  on  n  =  1, ...,  m  and  the  fact  that  -E{A}  =  /9';  the  third 
inequality  applies  the  condition  that  C  >  max(l,  (3')  and  m'^~^  >  m^~‘^  +  m  for  m  >  2.  This  concludes 
the  proof.  ■ 

Applying  Stirling’s  formula  and  Lemmas  1. 1,2,3,  we  can  obtain  the  conclusion  that  the  eigenvalue  distri¬ 
bution  of  SS^  in  the  shifted  model  is  identical  to  that  of  the  independent  model,  thus  assuring  the  assumption 
that  the  columns  of  S  can  be  regarded  as  independent  in  the  large  system  limit. 


Appendix  II 
Proof  of  Lemma  III. 2 

Proof:  From  the  definition  of  5a/,  we  have 

E{Ssif}  =  -  ^  (.E{S^5Sa}  -  E{5S^5Sa})  .  (24) 


We  consider  the  term  i?{5S^5Sa}  first.  It  is  easy  to  check  that  (recall  that  Ski  denotes  the  spreading  code 
of  user  k  along  path  /) 


1 

M 


1 

M 


M 

{m)E  {5bp5br} 

m=l 
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Jo,  ifp  ^  r 

1  if p  =  r  ’ 

where  p  =  \p\,  Q  =  mod  {i,L),  r  =  [f],  -s  =  mod(j,  L).  It  should  be  noted  we  applied  the  faet  that 
E  {5bp5br}  =  4Pe  in  the  seeond  equality. 

Aecording  to  Assumption  11.3,  the  spread  codes  are  mutually  independent  for  different  users  or  different 
paths.  Thus,  by  applying  the  strong  law  of  large  numbers,  we  have 


0,  if  (p,  q)  ^  (a  s) 

1,  if  (p,  q)  =  (a  s) 


Therefore,  we  have 


0,  if  i  #  j 

4A  if  1  =  7 

M  ’  ‘  J 


almost  surely,  as  M  ^  oo 


Similarly,  we  can  show  that 


J-.e{{S-6S)^ 


0,  if  i  ^  j 

7  _  j 

Ml  n  i  J 


almost  surely,  as  M  — >  cx) 


This  completes  the  proof. 


Appendix  III 


Proof  of  Prop.  III. 3 


Proof:  The  covariance  matrix  S/  is  given  by 


=  -^E{  Saa^fiS^S  }  “  {  S^(iSaa^(5S^(5S  } 

-  -^E  +  -^E  {(iS^(5Saa^(iS^(5S} 

-  E{Saf}E{Saf}^. 


The  elements  in  S^5Saa^5S^S  are  given  by 


M  M  KL  KL 


(S'^dSaa^bS'^S)..  ^EEEE^r  (p)(5sfc(p)sj(g)(5s;(g)afca;, 


p=l  7=1  k=l  1=1 
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where  Sj(p)  =  (p)S|-^1  njod(j,L)(p)’  namely  the  spreading  eode  (ineorporating  the  channel  symbol)  of 

the  mod(i,L)-th  path  of  user  \jr]  at  symbol  period  p,  6si(p)  =  (p)S|-^1  2,)(p)  and  a^  is  the  k- 

th  element  of  vector  a  and  equals  ap-|  To  compute  the  corresponding  expectation,  we  apply  the 

following  properties,  which  are  based  on  Assumption  II.4: 

•  Whenp  =  g,  if  I’^’l  =  \jr],  P(6sk{p)  ^  0,(5s/(g)  7^  0)  =  Pe,  since  Ssk{p)  and  6si{p)  are  determined 
by  the  same  decision  feedback; 

•  Whenp  =  g,  if  1"^]  7^  \jr],  P(6sk{p)  7^  0,(5s;(g)  7^  0)  =  P|,  since  6sk{p)  and  6si{p)  are  determined 
by  decision  feedback  from  different  users; 

•  When  p  ^  q,  P{6sk{p)  7^  0,(5s;(g)  7^  0)  =  Pg,  since  Ssk{p)  and  6si{p)  are  determined  by  decision 
feedback  from  different  symbol  periods; 

.  When  (5sfc(p)  7^  0,  6sk{p)  =  2sk{p). 

Thus  the  expectation  of  i  —  jth  element  of  S^5Sslsl^5S^S  is  given  by 

p|(S^(5Saa^(5S^S).J 

M  KL 

=  4Pe  EE  E  sj  {p)h{p)sj  {p)si{p)aka*i 

p=i  k=i  = 

M  KL 

+  ^Jip)^k{p)sj{p)si{p)aka*i 

p=i  k=i 

M  KL  KL 

+  4P2  ^  ^^§f(p)§fc(p)§J(g)si(g)afca; 

k=\  l=\ 

p,q=l 

p¥=<} 

=  T1+T2+T3, 

where  Ti,  T2  and  T3  represent  the  corresponding  three  summations,  respectively. 

Applying  the  strong  law  of  large  numbers  and  the  assumption  on  the  spreading  codes  that  {sj(p)}  are 
independent  for  different  values  of  i  or  p,  we  can  obtain  that,  as  M  — cx),  the  following  conclusions  hold 
almost  surely: 

4Pg  (^|a,|2  +  T  J2kX  k^i  lafcp)  ,  if  ^  =  J, 

^Pe  (1  +  jf)  a*a*,  if  i  7^  j  and  = 

0,  if  Til  ^  Til 
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4Pf  (l  +  a^a*, 


if  liWlil 


1 

M 


T2 


1 

JP 


n 


if  Til  =  Til 


Vi,i. 


We  can  apply  the  same  manipulation  and  obtain  that  E  |S^(5Saa^5S^(5S}  =  E  |(5S^(5Saa^(5S^S}  = 
\E  {6S'^6SsiSi^6S'^6S}  as  M  — 00.  Therefore,  we  can  obtain  (7)  since  the  sum  of  the  middle  three  terms 
in  (25)  is  zero  and  T3  cancels  E{6aLf}E{6aLf}^ .  ■ 

It  should  be  noted  that  the  above  analysis  is  also  valid  for  asynchronous  case  when  is  sufficiently  small. 
Similar  to  the  discussion  in  Section  III.  A,  we  can  remove  part  of  the  chips  in  the  first  and  the  last  symbol 
periods  to  obtain  a  similar  matrix  Sj^M-d^axxKL,  where  dmax  denotes  the  largest  time  offsets  of  different 
users,  measured  in  chips.  When  is  sufficiently  small  and  M  is  sufficiently  large,  we  can  ignore  the  terms 
scaled  by  P^  and  the  edge  effect  in  the  first  and  last  symbol  period.  Then,  we  have 

KL 

p{(S^<)Saa^<5S^S)^J  sfs,sfs^a,a;, 

^=1  \P=\P 

where  is  the  /c-th  column  of  matrix  S,  which  converges  to  Ti  as  M  — cx). 
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Fig  1.  CDMA  system  with  an  iterative  receiver 
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Variance 


Fig.  2.  Average  variance  of  channel  estimates  versus  the  coherence  time  M 
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Bit  Error  Rate 


Fig.  3.  Comparison  of  simulated  bit  error  rates  and  those  obtained  using  a  Gaussian  approximation 
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BER 


Fig.  4.  Performance  of  channel  codes  used  in  the  numerical  results,  where  the  input  SNR  =  ^ 
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max 


Fig.  5.  Maximum  load  of  systems  with  convolutional  codes 
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Fig.  6.  Maximum  load  of  systems  with  turbo  codes 
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Abstract 

A  game-theoretic  analysis  is  used  to  study  the  effects 
of  receiver  choice  on  the  energy  efficiency  of  multi-hop 
networks  in  which  the  nodes  communicate  using  Direct- 
Sequence  Code  Division  Multiple  Access  (DS-CDMA).  A 
Nash  equilibrium  of  the  game  in  which  the  network  nodes 
can  choose  their  receivers  as  well  as  their  transmit  powers 
to  maximize  the  total  number  of  bits  they  transmit  per  unit 
of  energy  is  derived.  The  energy  efficiencies  resulting  from 
the  use  of  different  linear  multiuser  receivers  in  this  context 
are  compared,  looking  at  both  the  non-cooperative  game 
and  the  Pareto  optimal  solution.  For  analytical  ease,  par¬ 
ticular  attention  is  paid  to  asymptotically  large  networks. 
Significant  gains  in  energy  efficiency  are  observed  when 
multiuser  receivers,  particularly  the  linear  minimum  mean- 
square  error  (MMSE)  receiver,  are  used  instead  of  conven¬ 
tional  matched  filter  receivers. 


1  Introduction 

In  a  wireless  multi-hop  network,  nodes  communicate 
by  passing  messages  for  one  another;  permitting  multi-hop 
communications,  rather  than  requiring  one-hop  communi¬ 
cations,  can  increase  network  capacity  and  allow  for  a  more 
ad  hoc  (and  thus  scalable)  system  (with  little  or  no  central¬ 
ized  control).  For  these  reasons,  and  because  of  their  poten¬ 
tial  for  commercial,  military,  and  civil  applications,  wireless 
multi-hop  networks  have  attracted  considerable  attention 
over  the  past  few  years.  In  these  networks,  energy  efficient 
communication  is  important  because  the  nodes  are  typi¬ 
cally  battery-powered  and  therefore  energy-limited.  Work 
on  energy-efficient  communication  in  these  multi-hop  net¬ 
works  has  often  focused  on  routing  protocols;  this  work  in¬ 
stead  looks  at  power  control  and  receiver  design  choices  that 
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can  be  implemented  independently  of  (and  thus  in  conjunc¬ 
tion  with)  the  routing  protocol. 

One  approach  that  has  been  very  successful  in  research¬ 
ing  energy  efficient  communications  in  both  cellular  and 
multi-hop  networks  is  the  game-theoretic  approach  de¬ 
scribed  in  [1,2].  Much  of  the  game-theoretic  research  in 
multi-hop  networks  has  focused  on  pricing  schemes  (e.g. 
[3, 4]).  In  this  work,  we  avoid  the  need  for  such  a  pricing 
scheme  by  using  instead  a  nodal  utility  function  to  capture 
the  energy  costs.  It  further  differs  from  previous  research 
by  considering  receiver  design,  as  [5]  does  for  cellular  net¬ 
works. 

We  propose  a  distributed  noncooperative  game  in  which 
the  nodes  can  choose  their  transmit  power  and  linear  re¬ 
ceiver  design  to  maximize  the  number  of  bits  that  they  can 
send  per  unit  of  power.  After  describing  the  network  and  in- 
ternodal  communications  in  Section  2,  we  derive  the  Nash 
equilibrium  for  this  game,  as  well  as  for  a  set  of  games  with 
set  receivers,  in  Section  3.  We  then  extend  the  asymptotic 
work  of  Tse  and  Hanly  [6]  to  fit  the  multihop  network  struc¬ 
ture  in  Section  4;  we  apply  this  in  Section  5  to  find  the 
Pareto  optimal  solution  in  an  asymptotically  large,  SINR- 
balanced  network.  Finally  we  present  some  numerical  re¬ 
sults  and  a  conclusion  in  Sections  6  and  7. 

2  System  Model 

Consider  a  wireless  multi-hop  network  with  K  nodes 
(users)  and  an  established  logical  topology,  where  a  se¬ 
quence  of  connected  link-nodes  I  e  L(k)  forms  a  route 
originating  from  a  source  k  (with  k  e  L{k)  by  definition). 
Let  m(k)  be  the  node  after  node  k  in  the  route  for  node  k. 
Assume  that  all  routes  that  go  through  a  node  k  continue 
through  m(k)  so  that  node  k  transmits  only  to  m(k).  Nodes 
communicate  with  each  other  using  DS-CDMA  with  pro¬ 
cessing  gain  N  (N  chips  per  bit). 

The  signal  received  at  a  node  m  (after  chip-matched  fil¬ 
tering)  sampled  at  the  chip  rate  over  one  symbol  duration 
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can  be  expressed  as 


k=\ 

where  p*.,  bj,,  and  are  the  transmit  power,  transmit¬ 
ted  symbol,  and  (binary)  spreading  sequence  for  node  h, 
hf^''  is  the  channel  gain  between  nodes  k  and  m;  and 
is  the  noise  vector  which  is  assumed  to  be  Gaus¬ 
sian  with  mean  0  and  covariance  cr^I.  (We  assume  here 
Pm  =  0.)  Assume  the  spreading  sequences  are  random, 
i.e.,  Sk  -  :^[fi  ■■■VnV,  where  the  v,’s  are  independent 
and  identically  distributed  (i.i.d.)  random  variables  taking 
values  {-l,-i-l)  with  equal  probabilities.  Denote  the  cross¬ 
correlations  between  spreading  sequences  as 

Pkj  —  ^k^j’  (2) 


noting  that  =  1  for  all  k. 

Let  us  represent  the  linear  receiver  at  the  mth  node  for 
the  kth  signature  sequence  by  a  coefficient  vector  cf'\  The 
output  of  this  receiver  can  be  written  as 

y  =  (3) 

=  ^  ^Jp-jhfbjCk^Sj  +  (4) 

j=kk 

The  signal-to-interference-plus-noise  ratio  (SINR),  of 
the  fcth  user  at  the  output  of  receiver  m(k)  is 

7k  =  - ^ - y-  (5) 


users;  and  /(•)  is  an  efficiency  function  that  closely  approx¬ 
imates  the  packet  success  rate.  This  efficiency  function  can 
be  any  increasing,  continuously  differentiable,  sigmoidal^ 
function  with  /(O)  =  0  and  /(-Hoo)  =  1.  See  [5]  for  more 
discussion  of  the  efficiency  function. 

Using  (7),  (6)  becomes 


^  ^u,f(yk) 
M 


(8) 


When  the  receiver  used  is  a  matched  filter  (MF)  (i.e. 
jjjg  received  SINR  is 

(rhk'^Sk  +  {sk'^Sjf 

Pkhf^" 

When  the  receiver  is  a  linear  minimum  mean-squared  er¬ 
ror  (MMSE)  receiver ,  the  filter  coefficients  and  the  received 
SINR  are  [7] 


(9) 

(10) 


and 


MMSE 


^/p'khf^ 


l+Pkhf^ 


(11) 


rf 


MSE 


:  Pkh 


k'^k. 


(12) 


where 


Each  user  has  a  utility  function  that  is  the  ratio  of  its 
effective  throughput  to  its  transmit  power,  i.e.. 


=  +  (13) 

j*k 


Uk  =  -.  (6) 

Pk 

Here,  the  throughput,  T^,  is  the  net  number  of  information 
bits  sent  by  k  (generated  by  k  or  any  node  whose  route  goes 
through  k)  and  received  without  error  at  the  intended  des¬ 
tination,  m{k),  per  unit  of  time.  (We  assume  that  all  the 
congestion  control  is  done  in  the  choice  of  routing.) 
Eollowing  the  discussion  in  [5],  we  will  use 

Tk  -  ^mjk)  (7) 

where  L  and  M  are  the  number  of  information  bits  and  the 
total  number  of  bits  in  a  packet,  respectively  (without  loss 
of  generality  assumed  here  to  be  the  same  for  all  users);  R  is 
the  transmission  rate,  which  is  the  ratio  of  the  bandwidth  to 
the  processing  gain  and  is  taken  for  now  to  be  equal  for  all 


When  the  receiver  is  a  decorrelator^  (DE)  (i.e.  C  = 
[ci  •  ■  ■  C/f]  =  S(S^S)“'  where  S  =  [si  ■  ■  ■  Sa:]),  the  received 
SINR  is 


I  m(k)^ 

Pkh^ 

O-^C^Ck  ' 


(14) 


Eor  any  linear  receiver  with  all  nodes’  coefficients  cho¬ 
sen  independently  of  their  transmit  powers  (including  the 
ME  and  DE),  as  well  as  for  the  MMSE  receiver. 


djk  _  n 
dpk  Pk ' 


(15) 


^  A  continuous  increasing  function  is  sigmoidal  if  there  is  a  point  above 
which  the  function  is  concave  and  below  which  the  function  is  convex. 
^Here,  we  must  assume  that  K  <  N. 
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3  The  Noncooperative  Power-Control  Game  As  long  as  the  users  all  have  the  same  efficiency  function, 


Let  0  -  I^TC,  {A^),  {Mj-lj  denote  the  noncooperative  game 
where  "TC  =  { 1 , . . . ,  /T)  and  At  -  [0,  Pmax]  x  is  the  strat¬ 
egy  set  for  the  feth  user.  Here,  P^ax  is  the  maximum  allowed 
power  for  transmission.  Each  strategy  in  At  can  be  written 
as  Ht  -  (Pk,  Ck),  where  pt  and  Ct  are  the  transmit  power  and 
the  receiver  filter  coefficients,  respectively,  of  user  k.  Then 
the  resulting  noncooperative  game  can  be  expressed  as  the 
maximization  problem  for  A:  =  \,. . .  ,K: 


LR 

max  ut  —  —  max 

at  M  Pk,<^k 


fijkiPk,  Ck)) 
Pk 


(16) 


where  yt  is  expressed  explicitly  as  a  function  of  pt  and  Ct- 
This  is  similar  to  the  noncooperative  power-control  game 
in  [5];  here,  however,  the  channel  gains  are  between  pairs 
of  nodes  rather  than  between  a  node  and  the  base-station. 

Since  the  choice  of  receiver  is  independent  of  the  trans¬ 
mit  power  and  /(■)  is  an  increasing  function,  the  analysis 
of  [5]  applies,  so  the  maximization  from  (16)  becomes: 


max 

Pt,Ct 


fiykipk,  Ck)) 

Pk 


/(maXcj7i(p^,Ci)) 

=  max - 

Pk  Pk 


(17) 


Note  that  the  MMSE  receiver  achieves  the  maximum  SINR 
amongst  all  linear  receivers,  so  that  if  a  Nash  equilibrium 
exists,  at  that  equilibrium  all  receivers  must  be  MMSE  re¬ 
ceivers.  Then  the  maximization  problem  becomes 


^(  MMSE(^^)) 

max - 

Pt  Pk 


(18) 


Let  0Q  —  [TC,  {[0, Pinax]), (M/t)]  denote  the  noncoopera¬ 
tive  game  that  differs  from  Q  in  that  users  cannot  choose 
their  linear  receivers  but  are  forced  to  use  the  receive  fil¬ 
ter  coefficients  [ci  ■  ■  ■  Ca:]  =  C  (which  may  be  a  function 
of  the  powers,  P).  The  resulting  noncooperative  game  can 
be  expressed  as  the  following  maximization  problem  for 
k^l,...,K\ 


,  ,  LR  firl^iPk))  ,,,, 

maxM^:  =  ma\ut(pk,Ct)  =  —  max -  (19) 

aj  Pt  M  Pt  Pt 

where  is  expressed  explicitly  as  a  function  of  pt-  Then 
the  maximization  problem  in  (18)  is  one  of  the  games  Qc 
when  C  is  chosen  to  be  the  MMSE  receivers. 

Eor  any  C  matrix  (or  C(P)  for  which  (15)  holds),  the 
utility  function  for  each  user  is  maximized  when 

Pk^m\n{Pt,pl]  (20) 

where  is  the  unique  positive  number  that  satisfies 

/(r:‘(pp)  =  7t(pl)f(Jt(pl)y  (21) 


r''(Pi)-...  =  r^"(py  =  r  (22) 


where  y*  is  the  unique  positive  number  that  satisfies 


f(y)  =  r/'(r  )■  (23) 

Einally,  since  is  quasi-concave^  in  pt,  we  can  use 
the  result  cited  in  [2,  Appendix  I]:  Qq  has  a  Nash  equilib¬ 
rium  and,  as  is  the  case  in  [5],  it  is  unique.  At  this  equilib¬ 
rium,  unless  there  is  a  node  k  with  >  Pt,  the  powers  are 
such  that  the  nodes  are  SINR-balanced  (i.e.  (22)  holds). 

Returning  to  the  game  0,  a  similar  result  holds:  there 
exists  a  unique  equilibrium  where  all  receivers  are  MMSE 
detectors  and,  if  the  power  limit  is  high  enough,  the  powers 
are  SINR-balanced. 

4  Asymptotically  Large  Systems:  Extend¬ 
ing  the  Tse-Hanly  Equations  to  Multi-Hop 
Networks 

Assume  that  the  channel  gains  are  independent.  That  is, 
in  the  asymptotic  regime  when  N,K  ^  oo  while  KjN  -  /3, 

(fflp 

the  interferers’  channel  gains,  for  all  m  +  k,m{k), 
are  iid  realizations  of  the  random  variable  G  with  pdf  fc, 
and  the  primary  channel  gains,  for  all  k,  are  iid 

realizations  of  the  random  variable  H  with  pdf  (where 
fuih)  =  OVh  <  0).  Let  q  -  ¥{m(f)  -  m(k)}  for  all  j  +  k. 

We  can  apply  results  from  [6]  to  analyze  the  nodes’ 
SINRs.  Then  we  find  a  probability  density  function  for  p 
such  that  in  an  asymptotically  large  system  where  all  nodes 
have  powers  distributed  by  this  function,  with  probability 
one  all  nodes  have  SINR  of  at  least  y  for  some  y.  If  this  dis¬ 
tribution  is  not  unique,  we  choose  the  one  that  minimizes 
the  nodes’  powers.  Eor  simplicity,  and  since  we  are  consid¬ 
ering  the  asymptotic  regime,  we  assume  that  the  distribution 
of  Pt  is  independent  of  all  channel  gains  except  for  . 

Eor  convenience  of  notation,  let 

Pt,^t 

nCC 

and  note  f{p,h)dp  -  fnih)  for  all  h.  Then  the  joint 
density  of  pt,  and  for  j  kis 

f  ,  (».w)2  (mjAp,  h,  g)  =  /p,h(p,  h)d(g  -  h)q 

+  fp,H(p,h)fc(8)(i-q)-  (24) 

Applying  the  results  from  [6],  when  the  receiver  at  node 
A:  is  a  matched  filter,  decorrelator,  or  MMSE  receiver,  the 
random  SINR  at  the  receiver  converges  in  probability  as 

function  is  quasi-concave  if  there  exists  a  point  below  which  the 
function  is  nondecreasing  and  above  which  the  function  is  non-increasing. 
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N,K  ^  oo  while  K/N  -  p.  These  asymptotic  SINRs  are  all  nodes  have  SINR  when  using  an  MMSE  receiver  of  at 
uniquely  described  by  the  equations  (where  j  +  k):  least  y  for  some  set  y.  Let  Q  =  inf {ph  :  f(p,  h)  >  0).  Then 


yMF  ^ 


CT^  +  I 


Pjli 


{m{k)Y 


a  <  1; 
a  >  1. 


(25) 


(26) 


and 


MSE 


Pkh' 


{m(k)y 


Sodgfp{p)fG{g)I{pg,  ,  yMMSE) 

(27) 


where /(a>,c)  = 

If  the  nodes  choose  their  transmit  powers  so  that  the 
SINRs  are  balanced,  the  following  theorem  determines 
what  SINRs  are  achievable  at  all  receivers  as  well  as  the 
minimum  transmit  powers  to  achieve  any  achievable  SINR 
when  the  nodes  use  the  MMSE  receiver,  under  the  assump¬ 
tions  listed  above. 


-><T^+ 13 
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Jo  Jo  Jo 
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dhfp^H(p,h)I(ph,  Q,y) 


XOO  /^OO  /^OO 

dgj^  dpj^  dhfp^H{p,h)fc{g)I{ph^,Q,7) 
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^00 

dp 
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dhfp,H(p,h)I(Q,  Q,y) 
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This  implies  that 
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1  H-  y 
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H  -tyG 
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h  +  yg 

(31) 


H  -HyG 


>  ycr^  >  0, 


(32) 


Theorem  4.1.  A  necessary  and  sufficient  condition  for  an 
SINR,  y,  to  be  achievable  is  for 


fyq 


1  H-  y 


+  Py(l  -q)E 


G 


H  +yG 


<  1. 


(28) 


When  (28)  holds,  each  user  can  achieve  the  desired  SINR, 
y,  and  the  minimum  power  solution  to  do  so  is  to  assign 
each  node,  k,  transmit  power 


Pk  -  Emmse  ,  yj 

1  ycr^ 

“  ■  1  -Pyq^^ -Py{\  -  q)E  ■ 


(29) 

(30) 


4.1  Proof  of  Theorem  4.1 


so/?y?i2_  +py(\-  q)E  [pj^]  <  1,  proving  necessity. 

When  (28)  holds,  it  is  easy  to  show  that  Pmmse  (h,  y)  is 
positive  for  all  primary  channel  gains,  h.  It  is  also  straight¬ 
forward  to  show  that  if  each  node,  k,  uses  transmit  power 

Pmf  ^^®®'  >  ^11  nodes  will  achieve  the  SINR  require¬ 

ment,  y,  finishing  the  proof  of  sufficiency. 

Einally,  consider  any  other  joint  distribution  of  powers 
and  primary  channel  gains  whose  marginal  distribution  for 
H  is  fu,  and  let  Q*  be  the  minimal  received  power  in  this 
distribution.  Then  by  exactly  the  same  argument  as  was 
used  in  the  proof  of  necessity. 


^  - Pyq^y  - PjO  - q)E[^] 
-  hPMMSE{h,y),  V/i  >  0. 


(33) 

(34) 


This  means  that  assigning  powers  according  to  Pmmse  does 
We  start  with  a  lemma  that  is  a  straightforward  conse-  j^deed  give  the  minimal  power  solution.  □ 

quence  of  the  definition  of  I{a,  b,  c). 


Lemma  4.2.  For  all  positive  real  numbers  ao,  a,  b,  c,  ao  <  a 
if  and  only  ifI(ao,  b,  c)  <  I(a,  b,  c). 

Then  the  proof  follows. 

Proof  To  show  necessity,  assume  that  there  is  a  pdf  /  with 
f(p,h)dp  -  fnih)  for  all  h,  such  that  in  an  asymptot¬ 
ically  large  system  where  all  nodes  have  powers  and  pri¬ 
mary  channel  gains  distributed  by  /,  with  probability  one 


5  A  Global  Optimization  Problem 


A  useful  global  optimization  problem  is 

K  r  K 


max 


Y,^kUk^-R  max^ 


k=l 


k=l 


akfijk) 

Pk 


(35) 


where  the  are  set  weighting  variables.  This  problem  is 
equivalent  to  finding  a  Pareto-optimal  solution  of  the  game. 
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According  to  [5],  even  in  the  special  case  of  a  cellular  sys¬ 
tem  where  L{k)  -  [k]  for  all  nodes  k  -  \,2,. . .  ,K  and  all 
nodes  are  transmitting  to  the  base-station,  “Pareto-optimal 
solutions  are,  in  general,  difficult  to  obtain.”  For  simplic¬ 
ity,  we  restrict  the  problem  by  requiring  that  the  solution 
is  “fair”:  all  nodes  have  equal  receiver  output  SINRs  (i.e. 
SINR-balancing),  so  7  =  yi  -  yi-  ■  ■  ■  Jk- 
With  this  assumption,  (35)  becomes 

^Rmax/(r)^^.  (36) 

M 

For  the  matched  filter,  we  can  apply  (5)  with  m  -  m(k) 
to  see  that  the  users’  SINRs  are  equal  if  and  only  if 


Any  jk  which  satisfies  ^  —  is  a  solution  to  (42). 

Then  the  power  for  user  k  to  achieve  the  SINR  y*  is 

pMMSE  ^  _ I _ _  ^43^ 

where  ((y)  is  the  mean  value  of  in  the  network.  Equal 
received  SINRs  amongst  the  users  is  achieved  with  mini¬ 
mum  power  consumption  when  =  K{y)  is  constant 

for  all  k  and 

2 

K(y)  =  - 7 — 7^ - r-  (44) 

1  -/3y{qj^  +(l-qX(y)) 


B -H  - -I- 1  m  p(7)  =  o-n 


where  B  is  a  Khy  K  matrix  with  entries  Bkj  —  p^j, 

D  is  K  hy  K  diagonal  matrix  with  diagonal  entries  0^.*.  = 
^(niW)  1  is  a  vector  of  K  ones. 

k  ’ 

The  SINR  that  maximizes  (36)  is  the  y  that  satisfies 


Then,  (36)  can  be  expressed  as 


o-^/(r)X^  (38) 

or[  Pkiy) 

oy  ^  Pk(y)  ^  pl(y)  dy 

(39) 

lere  Pkiy)  and  ^  {pk{y)\  the  kth  elements  of 


p(r)  =  (r^  B+  -  +  1  £>  1 


^  [p(r)]  ^  (P  {yB  +  {I  +  y)  D)  ^  D{yB  +  {\+y)D)  '  1. 
oy 


R 


fir) 

max-—, 
y  K(y) 


The  solution  to  maXy  ^  must  satisfy  =  0. 

Combining  this  with  (44)  gives  the  equation  that  must  be 
satisfied  by  the  solution  to  the  maximization  problem  in 
(45): 


fir)  =  r/'(r)  i  - 


-jS(l  -  q)yf(y) 


i-^-m-q)yf(r) 


If  ((y)  «  1,  then  the  equation  is  approximately  the 
same  as  in  the  cellular  case  [5]  with  K/N  — >  [iq.  Then, 
the  ability  to  use  multiple  hops  to  communicate,  and  there¬ 
fore  reduce  transmit  power,  has  similar  results  to  reducing 
the  system  load;  furthermore,  for  a  large  range  of  values  of 
Pq,  the  MMSE  target  SINRs  for  the  noncooperative  game 
and  for  the  Pareto-optimal  solution  are  close. 

6  Numerical  Results 


Eor  the  decorrelator,  it  is  easy  to  show  that  the  non- 
cooperative  results  are  equal  to  the  globally  optimal  results, 
since  the  users’  achieved  SINRs  are  independent  of  all  the 
powers  of  all  interferers. 

Einally,  for  the  MMSE  receiver,  we  can  apply  the  results 
from  Section  4.  In  a  large  system,  if  all  users  choose  their 
transmit  powers  based  on  the  values  of  for  m{j)  + 

m(k)  only  through  the  average  of  these  interference  gains 
and  if  we  use  the  assumptions  of  Section  4,  the  SINR  is 
approximated  by 

yvtMSE  ^  _ 

(42) 


Consider  a  multi-hop  network  with  /f  =  100  nodes  dis¬ 
tributed  randomly  in  a  square  500  meters  by  500  meters 
surrounding  an  access  point  in  the  center.  We  use  a  sim¬ 
ple  routing  scheme  where  all  nodes  transmit  to  the  closest 
node  that  is  closer  to  the  access  point  (or  the  access  point 
if  that  is  closest).  We  assume  that  each  packet  contains  100 
bits  of  data  and  no  overhead  (L  —  M  -  100);  the  trans¬ 
mission  rate  is  R  =  100  kb/s;  the  thermal  noise  power  is 
cr^  =  5  X  10“'®  Watts;  the  channel  gains  are  distributed 
with  a  Rayleigh  distribution  with  mean  0.3B“^,  where  d  is 
the  distance  between  the  transmitter  and  receiver;  and  the 
processing  gain  is  N.  We  use  the  same  efficiency  function 
as  [5],  namely  f{y)  =  (1  -  6“^)*^. 

Table  1  shows  the  average  utility  for  four  representative 
sets  of  randomly  chosen  spreading  sequences,  one  for  each 
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of  =  50,  100,  200,  and  300,  comparing  the  mean  utility 
under  the  various  power  choice  method  discussed  above. 
Table  2  shows  the  target  SINRs  for  the  socially  optimal  re¬ 
sults  displayed  in  Table  1 . 


MF 

DE 

MMSE 

A  =  50 
non-coop, 
soc.  opt. 

0 

2.025  X  10-'4 

1.198  X  10'° 

1.199  X  10'° 

N  =  100 
non-coop, 
soc.  opt. 

0 

1.050  X  lO-'* 

1.095  X  10" 
1.095  X  10" 

1.417  X  10'° 
1.417  X  10'° 

N  =  200 
non-coop, 
soc.  opt. 

0 

2.512  X  10-'° 

7.459  X  10° 
7.459  X  10° 

1.476  X  10'° 
1.476  X  10'° 

N  =  300 
non-coop, 
soc.  opt. 

0.2056 

1.351  X  10° 

1.001  X  10'° 
1.001  X  10'° 

1.493  X  10'° 
1.493  X  10'° 

Table  1.  Mean  utilities  for  four  representative 
sets  of  spreading  sequences. 


N 

MF 

DE 

MMSE 

50 

0.87 

6.39 

100 

1.31 

6.47 

6.43 

200 

0.99 

6.47 

6.45 

300 

5.03 

6.47 

6.46 

Table  2.  Socially  optimal  SINRs  for  the  same 
four  representative  sets  of  spreading  se¬ 
quences. 

The  socially  optimally  implemented  MF  receiver  per¬ 
forms  poorly  in  heavily-loaded  systems,  while  the  non- 
cooperative  implementation  fails  to  achieve  non-zero  utility 
except  in  the  case  with  the  lightest  load.  Even  in  the  case 
where  f5  -  1/3,  the  mean  utility  for  the  socially  optimal  MF 
receiver  is  less  than  a  tenth  of  the  MMSE  receiver’s  mean 
utility.  Using  the  DE  receiver  (for  which  we  require  that 
K  >  N),  as  was  noted  in  Section  5,  there  is  no  difference 
between  the  non-cooperative  and  socially  optimal  results: 
both  cases  have  the  same  target  SINR  and  thus  the  same 
mean  utility.  For  the  MMSE  receiver,  this  difference  be¬ 
tween  the  mean  utility  in  the  non-cooperative  and  socially 
optimal  implementations  is  very  small.  Finally,  the  DE  and 
MMSE  receivers  both  significantly  outperform  the  MF  re¬ 
ceiver  in  all  four  of  these  cases.  There  is,  however,  a  price  to 
pay  in  using  the  better-performing  receivers:  these  receivers 
require  more  information  at  every  node  as  well  as  signifi¬ 
cantly  more  computation.  These  issues  will  be  addressed 
further  in  later  research. 


7  Conclusion 

We  have  analyzed  the  cross-layer  issue  of  energy- 
efficient  communication  in  multi-hop  networks  using  a 
game  theoretic  method.  Focusing  on  linear  receivers,  we 
have  derived  the  transmit  power  levels  that  results  in  a  Nash 
equilibrium  for  multiple  receiver  designs,  showing  that  at 
this  equilibrium  the  users  are  SINR-balanced.  We  then  gen¬ 
eralized  the  important  asymptotic  work  of  Tse  and  Hanly 
to  allow  for  the  case  where  users  and  their  interferers  may 
be  transmitting  to  different  locations,  keeping  the  cellular 
example  as  a  special  case.  We  applied  these  asymptotic  re¬ 
sults,  as  well  as  exact  results  for  the  MF  and  DE  receivers, 
to  find  the  equations  for  the  SINR-balanced  Pareto-optimal 
solution.  We  showed  that  the  MMSE  receiver  is  the  optimal 
receiver  and  that  in  many  cases  the  non-cooperative  MMSE 
receiver  results  are  quite  close  to  the  socially  optimal  re¬ 
sults. 
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